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1.0  Objectives  and  their  evolution 

Light  atom  doped  solid  hydrogen  is  one 
of  the  most  promising  candidates  for  a 
practical  HEDM  that  can  be  used  to  improve 
the  performance  of  both  civilian  and  military 
laimch  vehicles.  An  important  goal  of  research 
in  this  field  is  to  determine  the  maximum 
concentrations  that  can  be  achieved  while 
maintaining  a  high  energy  density  of  the  tliel. 
Answering  this  question  requires  improved 
imderstanding  of  hydrogen-dopant  and  dopant- 
dopant  interactions  and  the  mobility  of  dopants 
in  both  liquid  and  solid  hydrogen.  Our 
contribution  to  this  effort  has  been  the 
development  and  exploitation  of  methods  to 
produce  and  spectroscopically  characterize 
metal  doped  quantum  clusters  made  of 
hydrogen  and/or  helium.  Our  work 
complements  ongoing  research  on  doped  solid 
hydrogen  at  Edwards  AFB. 

Using  spectroscopy  as  a  diagnostic  tool 
our  program  aimed  at  the  study  of  metal-doped 
H2  clusters  in  order  to  investigate  the 
possibility  that  certain  light  atoms  or  their 
dimers  may  have  small,  but  significant, 
barriers  to  Ae  formation  of  covalently  bound 
larger  clusters.  At  the  temperature  of  solid 
hydrogen  fuels,  the  presence  of  even  small 
activation  barriers  (on  the  order  of 
0.2  kcal/mole)  will  reduce  the  rate  of  chemical 
bond  formation.  As  the  rate  of  bond  formation 
becomes  lower,  higher  dopant  concentrations 
could  be  realized  in  condensed  hydrogen  fuel. 

The  first  approach  we  have  explored  is 
to  exploit  small  barriers  that  may  arise  in  the 
recombination  of  some  atoms  due  to  the  need 
to  decouple  spin-orbit  interactions  in  order  to 
have  atomic  orbitals  that  can  participate  in 
bonding.  This  approach  arose  firom  our 
previous  work  on  the  dynamics  of  excited 
states  of  alkali  atoms  on  the  surface  of  helium 
nanodroplets,  where  such  barriers  were  shown 
to  strongly  influence  the  dynamics.  A  second 
approach  is  to  build  up  van  der  Waals  clusters 
of  alkaline  earth  atoms,  such  as  Be  or  Mg, 
which  are  closed  shell  species  and  thus  may 
have  small  barriers  to  ‘metallization.’  By  thus 


keeping  the  size  of  the  cluster  below  the 
metallization  limit,  the  energy  contents  of  the 
fuel  is  optimized.  A  third  approach  is  to 
saturate  the  interactions,  thus  kinetically 
preventing  the  formation  of  covalent  bonds. 

A  recent  theoretical  investigation  of 
boron-doped  H2  matrices  by  Voth,  Alexander 
and  collaborators  has  shown  that  at  a 
concentration  of  a  few  percent  a  boron-doped 
H2  matrix  should  be  stable.  On  the  other  hand 
experimental  work  by  Fajardo  and 
collaborators  at  the  AFRL  (Edwards  AFB)  has 
failed  to  show  that  in  the  case  of  A1  doping 
such  concentration  can  be  reached  as 
catastrophic  atom-atom  recombination  seems 
to  set  in  already  at  concentrations  below  the 
1%  level.  However,  it  is  quite  likely  that  the 
diffusional  processes  that  lead  to  the 
imstablization  of  the  matrix  occur  at  the 
moment  the  matrix  is  formed.  This  is  because 
the  high  H2  flux  needed  to  obtain  transparent 
matrices  (i.e.,  matrices  that  can  be  probed) 
leads  necessarily  to  “warm”  surface  conditions 
which,  in  turn,  allow  for  dopant  diffusion. 
(The  latter  is  also  favored  by  the  fact  that,  at 
present,  the  A1  atom  flux  comes  from  a  hot 
metal  source). 

The  study  of  metal-doped  H2  clusters 
(with  or  without  a  liquid  He  coating)  in 
progress  in  our  laboratory  aims  a  better 
imderstanding  of  the  location  of  the  dopant  in 
the  cluster  towards  an  evaluation  of  the 
feasibility  to  produce  higher  (~5%)  metal- 
doped  matrices  by  cluster  deposition  instead  of 
beam  (gas)  deposition  as  used  at  present  at  the 
AFR  Laboratory. 

Since  the  results  obtained  in  the  early 
part  of  the  reporting  period,  on  the  visible/uv 
spectroscopy  of  Mg  and  A1  atom-doped  He 
clusters  had  shown  the  limits  of  this  type  of 
diagnostic  (too  broad  lines  instead  of  strong 
inelastic  probe  atom/matrix  interactions)  we 
have  concentrated  in  the  second  half  of  the 
period  on  obtaining  IR  spectra  of  molecule 
doped  H2  clusters  and  on  interpreting  the 
results  obtained  last  year.  We  have  made 
strong  progress  on  the  interpretation  side  so 
that  a.)  we  can  now  predict  with  confidence  the 
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location  of  a  guest  atom  in  a  quantum  cluster 
made  of  He  or  H2  and  b.)  we  understand  both 
the  multiplicity  and  the  lifetimes  of  the  spectra 
measured  last  year.  We  have  instead  failed  in 
obtaining  IR  spectra  of  molecules  embedded  in 
H2  clusters.  This  is  likely  due  to  the  fact  that 
we  are  limited  by  the  5000-6700  cm'' 
frequency  range  of  our  tunable  laser  system 
which  gives  us  easy  access  to  the  overtone 
region  of  CH  chromophores.  It  is  likely  that 
these  overtone  excitations  relax  too  rapidly  in 
H2  clusters  producing  spectral  features  too 
broad  to  be  detected  within  our  ultrahigh- 
resolution  apparatus. 

Two  years  ago  we  have  been  fortunate 
to  obtain  a  DURIP  grant  with  which  we  have 
purchased  a  new  laser  system  (an  Indigo  laser 
manufactured  by  Positive  Light)  which,  in 
combination  with  a  laser  built  in  our 
laboratory,  was  designed  to  obtain  high 
resolution  Raman  spectra  of  both  metal-doped 
and  undoped  clusters.  Unfortunately  when  the 
Indigo  laser  arrived  about  14  months  ago  it  did 
not  meet  our  expectations  with  respect  to 
linewidth  and  stability.  On  the  other  hand,  as  it 
was  satisfactory  with  respect  to  the  repetition 
rate  and  the  average  power  levels  that  it  could 
reach,  and  because  there  was  not  an  alternative 
solution  available  on  the  market,  we  decided  to 
tough  it  out  and  fix  the  laser  ourselves.  This 
effort,  which  resulted  to  be  somewhat  more 
difficult  than  expected,  lasted  approximately 
one  year;  but  now  (see  also  the  DURIP  grant  # 
F49620-00-1-0192  final  report  filed  at  the 
same  time  as  the  present  one)  we  are  happy  to 
report  that  the  laser  delivers  an  average  power 
of  1 .5  Watts  at  the  repetition  rate  of  1  kHz  with 
pulse  length  10  ns  and  bandwidth  of  725- 
830  ran.  This  level  of  performance,  which 
makes  our  system  a  rather  unique  piece  of 
equipment,  should,  barring  unforeseen 
circumstances,  allow  us  in  the  next  few  months 
to  obtain  Raman  spectra  of  both  doped  and 
undoped  hydrogen  cluster  beams.  Raman 
excitation  of  H2  molecules  should  show  much 
narrower  spectral  features  which,  by  shifts  and 
splittings,  should  also  reveal  the  presence  and 
location  of  eventual  metal  atom  dopants.  In 


this  way  the  problems  described  above  should 
be  eliminated  and  a  rich  selection  of 
experimental  possibilities  should  become 
available. 


2.  Accomplishments/New  findings 

Contents: 

2.1_  Time  resolved  studies  of  complexes 
involving  alkali  atoms. 

2.2_  Solvation  of  Mg  atoms  in  superfluid 
helium  droplets. 

2.3_  Solvation  of  A1  atoms  in  superfluid 
helium  droplets. 

2.4_  Theory  of  intermolecular  forces  and  its 
application  to  light-metal  atom 
complexes. 

2.5_  Progress  towards  assembling  a  facility 
for  stimulated  Raman  spectroscopy  in 
doped  and  undoped  hydrogen  clusters. 

2.1  Time-resolved  studies  of  complexes 
involvine  alkali  atoms 

2.1.1  Dynamics  of  alkali-helium  excimers 

The  formation  of  a  Na^-He  excimer 
upon  optical  excitation  of  a  Na  atom  on  the 
surface  of  a  helium  droplet  has  been  previously 
demonstrated  in  our  laboratory.’’^  The 
emission  is  red-shifted  from  the  Na  atomic 
lines  and  can  be  quantitatively  modeled  from 
the  potential  curves  of  J.  Pascale.^  Using  time- 
correlated  single  photon  counting  we  have 
measured  the  onset  of  fluorescence  from  the 
Na*-He  excimer;  depending  on  which  Na 
excited  state  we  access,  we  see  different  rates 
of  fluorescence.  The  rise  time  of  excimer 
fluorescence  after  excitation  of  the  Na  3Pi/2 
line  was  fovmd  to  be  700  ps,  while  after 
exciting  the  Na  3P3/2  line  the  rise  time  was 
found  to  be  70  ps. 

The  difference  in  rise  time  has  been 
attributed  to  spin-orbit  mixing  of  the  ^ni/2 
surface  of  NaHe  with  the  surface, 
generating  a  small  (0.7  cm'')  barrier  on  the 
^ni/2  (which  asymptotically  corresponds  to  the 
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Table  1.  Rise  and  fall  times  o1 

'  intersystem  crossing 

products  due  to  quartet  Na  trimer  excitation.  The 
columns  v  and  j  indicate  the  quantum  numbers  of  the 
excited  trimer  bands  (listed  in  order  of  descending 
energy). 

Na  J=l/2  level)  but  no  barrier  on  the  ^113/2 
curve  (asymptotically,  the  J=3/2  level).  Even 
though  tliis  barrier  is  small,  there  will  be 
significant  effects  on  the  dynamics  due  to  the 
extremely  cold  He  droplet  environment.  A 
similar  spin-orbit  induced  barrier  as  been 
observed  for  Ag*-He  (in  bulk  liquid  helium).^ 

We  have  also  characterized  the 
formation  of  the  K*-He  excimer  formed  after 
optical  excitation  of  K  atoms  on  the  surface  of 
a  He  droplet.  Again,  the  dispersed  emission  is 
red-shifti  and  can  be  modeled  using  tiie  KHe 
surfaces  of  Pascale.^  The  fit  is  good  except  in 
the  red  tail  which  may  be  due  to  excessive 
stiffiiess  in  the  upper  state  repulsive  wall,  or  to 
possible  higher-order  K*-Hen  exciplexes. 

In  order  to  compare  to  Na*-He,  we 
have  measured  the  onset  of  the  red-shifted 
emission  from  excitation  of  K  atoms.  We 
record  a  rise  time  of  the  excimer  fluorescence 
of  200  ps,  quite  smaller  than  for  the 
comparable  Na*-He  case.  One  would  expect 
an  increase  in  the  size  of  the  barrier  on  the 
surface  (because  of  the  larger  spin-orbit 
splitting),  leading  to  a  longer  time  before  onset 
of  fluorescence.  However,  the  spin-orbit 
splitting  is  not  the  only  factor.  One  must  also 
consider  the  long  range  dispersion  interaction 
to  properly  predict  the  presence  and  magnitude 
of  a  barrier  along  the  spin-orbit  decoupled 
potential  energy  surface. 


V 

j 

Fall  Time 
(±0.1  ns) 
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H9i 
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Table  2.  The  lifetimes  of  the  v,  j  vibrational  levels  of 
the  2^E’  state  of  Nas. 


In  the  case  of  KHe,  the  dispersive 
forces  are  much  stronger  than  for  NaHe.  For 
KHe,  the  Ce  coefficient  is  38.9  atomic  units, 
but  only  25.5  a.u.  for  NaHe^.  Any  potential 
barrier  that  is  formed  through  spin-orbit 
interaction  is  tempered  by  the  relative  strength 
of  the  long  range  attractive  forces,  resulting  in 
a  reduced  barrier  and  faster  onset  of  excimer 
fluorescence.  Two  papers  have  been  published 
on  this  subject  in  tl;e  Journal  of  Chemical 
Physics^’’. 

2.1.2  Dynamics  of  alkali  quartet  trimers 

We  have  examined  the  curve  crossing 
undergone  by  spin-polarized  Naa  upon  optical 
excitation  to  various  vibronic  bands  of  the  2‘*E’ 
electronic  state.*’  The  excited  quartet  state 
crosses  to  a  predissociative  doublet  state, 
opening  two  emission  chaimels  from  the 
doublet  trimer:  dissociation  into  a  singlet 
dimer  and  an  excited  atom,  or  an  excited  dimer 
and  a  ground-state  atom.  In  both  channels, 
fluorescence  is  blue-shifted  from  the  excitation 
frequency.  K3  goes  through  a  similar  spin-flip 
reaction  when  electronically  excited  in  the 
quartet  manifold. 

For  both  species,  we  have  measured  the 
intersystem  crossing  times  as  a  function  of 
excited  quartet  vibronic  band.  The  crossing 
time  is  obtained  by  timing  the  onset  of 
fluorescence,  using  filters  to  select  the  desired 
channel.  The  rise  times,  after  excitation  of 
different  vibronic  bands,  for  both  the  Naa 
singlet  (B->X)  and  Na  atomic  (3P->3S) 
emission  are  listed  in  Table  1,  along  with  the 
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measured  lifetimes  of  the  (B)  I'flu  state  and 
the  atomic  3P  state.  The  lifetimes  of  both 
states  correspond  well  with  the  gas-phase 
literature  values,'®  indicating  the  products  of 
the  spin-flip  are  no  longer  attached  to  the 
droplet  surface. 

The  onset  time  represents  the  time 
needed  to  populate  the  two  emission  channels, 
which  is  also  the  time  needed  for  intersystem 
crossing  from  the  quartet  to  the  doublet 
manifold.  Since  the  branching  ratio  between 
the  two  channels  is  constant  with  respect  to  the 
quartet  vibronic  band  excited,  we  can  conclude 
diat  there  is  a  single  access  window  into  the 
doublet  manifold,  located  near  16030  cm  '.  A 
paper  has  been  published  in  the  Journal  of 
Chemical  Physics  on  this  subject". 

The  resonant  fluorescence  of  quartet 
Na3  has  also  been  measured  in  time.  This  is 
the  excitation  to  the  2'‘E’  state,  which  then 
fluoresces  back  to  the  ground  quartet  state, 
I'^Aa’.  All  the  rise  times  are  under  40  ps, 
which  is  the  limit  of  our  instrument,  while  the 
fall  times  are  several  hvmdred  nanoseconds  (see 
Table  2)  and  biexponential.  As  excitation 
energy  is  increased,  the  decay  time  of  the 
resonant  fluorescence  decreases.  This  matches 
the  trend  of  the  rise  time  of  both  the  atomic  and 
dimer  fluorescence.  Further,  that  the  lifetimes 
of  the  excited  quartet  states  closely  match  the 
onset  time  of  the  products  provides  evidence 
for  the  minimal  influence  of  the  helium 
droplet,  since  the  majority  of  depopulation 
occurs  through  intersystem  crossing  and 
radiative  decay.  From  the  resonant 
fluorescence  spectra  we  have  obtained 
substantial  information  on  three-body  forces  in 
low  Ionization  Potential  atoms  and  published  a 
paper  on  this  subject  in  the  Journal  of 
Chemical  Physics'^ 

2.1.3  Dynamics  of  triolet  potassium  dimers 

Along  with  curve-crossing  in  alkali 
timers,  we  have  also  observed  in  the 
l^rig^-l^Xu^  transition  of  K2  not  only  atomic 
4P— >4S  emission  but  also  K2  B('nu)— >  and 
A('Zu^->  X('Eg^  dimer  emission.  This 


emission  is  again  blue-shifted  and  indicative  of 
an  intersystem  crossing  mediated  by  the 
environment  of  the  helium  droplet.  While 
scanning  over  the  1  ng<— 1  Zu  transition,  we 
selectively  collected  the  B->X  and  A->X 
singlet  emissions.  The  fluorescence  from 
molecular  emission  is  shifted  83  cm’'  to  the  red 
of  the  fluorescence  from  atomic  emission  after 
the  same  excitation. 

Time-resolved  studies  were  made  of  the 
fluorescence  arising  from  the  l^Ilg^— l^Zu^ 
excitation,  particularly  the  rise  times  of  the 
emission  from  all  three  species  produced.  For 
both  singlet  states,  the  fluorescence  starts  65  ps 
after  excitation,  while  the  atomic  fluorescence 
occurs  in  less  than  50  ps. 

The  potential  surface  of  the  l^ITg  state 
of  K2  contains  four  vibrational  levels  held 
behind  a  small  barrier.  There  is  an  83  cm*' 
separation  between  the  first  and  fourth  levels. 
The  lifetime  of  these  levels  decreases  as  the 
vibrational  number  increases.  We  can 
conclude  that  excitation  into  the  higher 
vibrational  levels  of  the  triplet  excited  state 
causes  dissociation  and  atomic  fluorescence. 
And  excitation  into  the  lower  vibrational  levels 
(primarily  v’=0,  the  only  level  with  a  lifetime 
>60ps)  will  lead  to  an  intersystem  crossing 
and  singlet  dimer  fluorescence.  On  this  subject 
we  have  published  a  paper  in  the  Discussions 
of  the  Faraday  Society". 

2.2  Solvation  of  Mg  atoms  in  superfluid 
helium  droplets  < 

2.2.1  Mg  atom  excitation  spectra 

Our  measurement  of  the  3'Pi*^  4-  3'So 
transition  of  Mg  atoms  in  helium  nanodroplets 
shows  evidence  that  this  atom  is  solvated  in  the 
interior  of  the  droplets  (see  Figure  1).  The 
absorption  is  broadened  by  the  droplet 
(FWHM  of  660  cm’')  and  blue-shifted  from  its 
gas-phase  transition  (Ay  =  307  cm"'  and 
642  cm'').  The  spectrum  shows  a  splitting  into 
two  components  that  are  attributable  to 
quadrupolar  deformations  of  the  helium  cavity 
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Figure  1.  (upper  panel)  The  3’Pi®  <-  3'So  transition  of 
Mg  atoms  solvated  in  helium  nanodroplets.  The  best  fit 
to  the  data  by  two  Gaussian  lineshapes  is  also  shown 
(solid  lines).  Also,  the  same  transition  in  bulk  liquid 
helium  (dots).  Qower  panel)  The  same  transition  in 
Ca;  both  in  helium  nanodroplets  (solid),  and  in  bulk 
liquid  helium  (dots). 

surrounding  the  atom.  This  splitting  has  been 
seen  by  Kinoshita  et  al.  in  Da 
(n^  P3/2  <-  n^  S1/2)  excitation  of  Cs  (n  =  6)  and 
Rb  (n  =  5)  atoms  in  bulk  liquid  helium,  as  well 
as  our  studies  of  Al‘^  atoms  solvated  in  liquid 
helium  droplets. 

As  the  droplet  spectrum  matches  the 
spectra  taken  in  a  bulk  liquid  helium  matrix'*, 
the  similarity  of  the  two  spectra  speaks  in  favor 
of  solvation.  In  contrast,  alkali  atoms  (on  the 
helium  droplet  surface)  exhibit  only  small 
shifts  from  their  gas-phase  frequencies.  By 
comparing  Ca  on  a  droplet'^  with  Ca  in  bulk 
liquid  helium'*,  we  can  see  that  it  holds  an 
intermediate  position  between  solvated  and 
surface  locations,  having  only  one-half  of  the 


broadening  and  one-third  of  the  shift  as  the 
bulk  helium  spectra. 

2.2.2  Me  atom  time-resolved  measurements 

Along  with  the  excitation  spectra,  we 
have  temporally  dispersed  the  fluorescence 
emitted  after  S'?!®  3' So  excitation  of  the 

solvated  Mg  atoms  (see  Figure  2).  Both  the 
gas-phase  and  the  droplet-solvated  Mg  spectra 
are  well  fit  by  a  model  employing  a  single 
exponential  rise  and  a  single  exponential  fall. 
The  fall  time  (i.e.,  lifetime)  of  2.05  ±0.05  ns 
found  for  the  3'Pi°  state  emission  of  gas-phase 
Mg  agrees  well  with  the  literature  value 
(1.99  ±  0.08  ns)'^  while  the  rise  time  fit  to  the 
data  is  instrument-limited  (i.e.,  <  50  ps).  For 
solvated  Mg  atoms  in  the  droplet,  the  lifetime 
is  found  to  be  2.39  ±  0.05  ns,  which  is  20% 
longer  than  the  gas-phase  value. 

The  lengthening  of  the  lifetime  for  the 
solvated  atom  can  be  ascribed  to  two  sources, 
the  first  being  the  red-shift  of  the  emission 
frequency.  Since  the  rate  of  spontaneous 
emission  scales  with  the  cube  of  the  frequency, 
a  lower  emission  frequency  will  give  a  longer 
lifetime.  Assuming  a  red-shift  similar  to  that 
of  Mg  atoms  solvated  in  bulk  liquid  helium'^, 
which  is  3.2±0.2nm,  accounts  for  a  3.5% 
increase  in  the  emission  lifetime. 

The  second  cause  of  the  increased 
lifetime  is  the  helium  solvent  itself.  The 
dielectric  response  of  the  solvent  will  result  in 
an  induced  transition  dipole  that  will  interfere 
with  the  transition  dipole  of  the  chromophore. 
By  summing  over  all  the  dipoles  induced  by 
the  radiating  dipole,  we  can  write  the  net 
induced  dipole  as  the  following  integral: 

=//-2a-|pf)P2(cos0)-r'Vv 

where  r  is  the  distance,  a  is  the  polarizability 
volume  of  He  (0.206  A^),  and  d  is  the  angle 
between  the  chromophore  transition  dipole  and 
the  displacement  to  a  generic  point  in  the  He 
cluster.  Taking  into  account  the  anisotropic  He 
density  around  the  excited  Mg  atom,  shown  in 
Figure  3,  (with  high  density  at  the  node  of  the 
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Figure  2.  Log  plot  of  the  time  evolution  of  gas-phase 
(filled  squares)  and  droplet-solvated  (open  circles)  Mg 
atoms  upon  3'Pi®4-3'So  excitation.  The  droplet- 
solvated  data  are  normalized  to  the  gas-phase  data  for 
visual  comparison.  The  solid  lines  are  the  best  fits  to  the 
data. 

p  orbital,  which  is  parallel  to  the  transition 
moment),  the  induced  moment  is  6.2%  of  the 
bare  Mg  transition  moment.  This  leads  to  a 
predicted  lifetime  lengthening  of  14%.  When 
combined  with  the  3.5%  from  the  emission 
frequency,  we  predict  a  lifetime  17.5%  longer 
than  the  gas-phase  value.  Our  experimental 
value  of  20  ±3%  adds  further  evidence  of 
solvation  by  the  helium  droplet.  These  results 
have  been  published  in  ref  20. 


3.3  Mg-He  potential  energy  surfaces 

Ancilotto  et  have  used  the  dopant- 
He  pair  potential  to  devise  a  relatively  simple 
model  that  treats  the  question  of  solvation. 
They  propose  that  solvation  can  be  predicted 
by  means  of  a  dimensionless  parameter  X, 
which  compares  the  gain  in  energy  due  to  the 
dopant-helium  interaction  against  the  energetic 
cost  of  creating  a  cavity  within  the  liquid. 
Formally, 

where  p  is  the  helium  density,  a  is  the  surface 
tension  of  the  liquid,  and  £  and  r^in  are  the  well 
depth  and  equilibrium  bond  distance  of  the 
dopant-He  pair  potential,  respectively. 
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Figure  3.  He  density  contour  plot  calculated  via 
Density  Functional  Theory  in  which  a  ground  state  Mg 
atom  is  constrained  to  be  in  the  center  {i.e.,  at  0,  0)  of  a 


300-atom 

lelium  droplet. 

- 0 - 

min(A) 

X 

IP(eV) 

Cs 

1.2 

7.7 

0.66 

3.9 

Rb 

1.4 

7.3 

0.73 

4.2 

K 

1.0 

7.7 

0.59 

4.3 

Na 

1.7 

6.5 

0.85 

5.1 

Li 

2.2 

6.1 

1.01 

5.4 

Ca 

3.7 

6.4 

1.8 

6.1 

Mg 

5.6 

5.3 

2.2 

7.6 

Table  3.  A,  values  fi'om  HFD  dopant-He  potentials. 
Threshold  of  solvation  should  be  between  Ca  and  Mg. 


Using  helium  density  functional 
calculations,  Ancilotto  et  al.  predict  that 
atoms  with  values  of  X,  <  1 .9  reside  on  the 
surface,  while  dopants  with  X  >  1.9  remain 
solvated.  The  model  successfully  predicts 
exterior  locations  for  alkali  atoms  (X  =  0.7)  and 
interior  locations  for  Ag  (Xs5)  and  SFs 
(Xsl9).  Table  3  lists  the  parameters  for 
several  alkali  and  alkaline  earth  solvated 
atoms. 

We  have  evaluated  the  potential  energy 
surfaces  of  systems  at  the  border  region  of  the 
Ancilotto  model.  In  so  doing,  we  have 
produced  the  (X)1^E  surfaces  of  MgHe  and 
CaHe  using  the  Hartree-Fock  Damped 
Dispersion  ansatz,  a  widely  used  semi- 
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Figure  4.  Potential  surfaces  of  MgHe. 

empirical  method  which  has  been  found  to  be 
quite  accurate  for  atom  combinations  involving 
atoms  of  the  first  and  second  rows  of  the 
periodic  table““. 

Figure  4  illustrates  the  MgHe 
(X)l^Z  state  as  calculated  by  several  methods. 
Well  depth  e  and  equilibrium  intemuclear 
distance  rmin  for  these  surfaces  are  reported  in 
Table  4.  Of  the  surfaces  available  for  both 
dopants  only  the  HFD  and  CCSD(T)  yield  X 
values  for  Mg  and  Ca  that  follow  the  expected 
trend.  The  CCSD(T)/Ames  MgHe  potential 
gives  a  X,  value  of  2.6,  correctly  predicting 
solvation;  for  Ca,  a  X  value  of  2.2.  This  is  just 
over  the  threshold  for  solvation  predicted  by 
the  Ancilotto  model,  establishing  the 
“borderline”  nature  of  Ca  attached  to  liquid 
helium  droplets.  One  way  to  picture  the 
“borderline”  case  is  a  half-submerged  atom 
with  ripples  of  helium  density  passing  over  the 
top. 

Helium  density  fimctional  studies  have 
been  used  as  well  to  explore  this  situation.  Our 
density  functional  studies  were  done  using 
code  for  helium  with  impurity  atoms  in  which 
the  dopant  atom  is  treated  as  a  classical  object 
(/.e.,  devoid  of  zero  point  energy)  in  a  quantum 
liquid^*’^’.  Previous  calculations  employing 
these  programs  have  been  found  to  be  in  good 
agreement  with  Monte  Carlo  simulations^®. 

In  these  studies,  the  solvation  energy 
Esoivated  of  the  guest  atom  inside  a  Hen  droplet 


e(cm‘) 

X 

Mg 

Ca 

Mg 

Ca 

Mg 

Ca 

SCF/CI 

2.3 

2.9 

5.6 

7.4 

1.4 

2.1 

A-T-T 

7.7 

10.3 

4.7 

5.1 

OO 

5.5 

HFD 

3.9 

2.6 

5.3 

6.4 

2.22 

1.82 

CCSD(T)/Ames 

4.8 

3.3 

5.1 

6.0 

2.65 

2.16 

CCSD(T)/UD 

4.0 

5.2 

2.2 

SAPT 

5.8 

5.0 

2.7 

Table  4.  Comparison  of  the  well  depth  (e),  intemuclear 
distance  (t^),  and  the  X  paramter  of  Ancilotto  et  al^' 
for  the  l^X  surfaces  of  MgHe  and  CaHe  as  calculated  by 


the  SCF/CP,  A-T-T^®,  HFD“,  CCSD(T)/UD^*,  SAPT^*, 
CCSD(T)/Ames^’  methods.  _ 


Dopant 

Esurface(K) 

Esolvated(Kj 

Na 

-11 

+46 

Ca 

-7.5 

+80 

Mg-HFD 

-21 

-0.2 

Mg-CCSD(T)/UD 

-16 

+2.7 

Mg-SAPT 

-32 

-41 

Table  5.  Energies  from  density  functional  calculations 
for  various  dopants  in  and  on  a  helium  droplet  of  300 
helium  atoms.  All  energies  are  normalized  to  the 
droplet-dopant  asymptote. 

(n  =  300)  is  compared  with  the  energy  of  the 
same  atom  resident  on  the  surface  of  a  droplet 
of  the  same  size.  The  dopant  atom  is  treated  as 
an  external  potential  around  which  the  helium 
density  is  allowed  to  adjust  based  upon 
minimization  of  the  total  energy  of  the  system. 
Esurface  is  calculated  with  the  dopant  fixed  at  16 
A  fi-om  the  center,  corresponding  to  the  radius 
of  the  300-atom  droplet,  while  Esoivated  has  the 
dopant  constrained  to  be  at  the  center  of  the 
droplet. 

The  DFT  results  for  Na  and  Ca,  in 
Table  5,  are  correctly  predicted  to  be  imstable 
when  solvated,  but  the  Mg  results  are 
unfortunately  not  right.  All  three  MgHe 
potentials  lead  to  minima  when  on  the  surface; 
only  the  SAPT  potential  (which  is  too  deep) 
also  has  a  minimum  in  the  interior.  This  is 
more  a  function  of  the  density  functional 
putting  too  much  emphasis  on  the  repulsive 
wall  than  of  the  SAPT  potential  being  correct. 
Looking  at  the  He  density  profiles,  it  is  clear 
that  the  SAPT  potential  creates  a  smaller  cavity 
in  the  helium  and  is  therefore  less  likely  to  be 
expelled  to  the  surface  (see  fig.  5). 
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Figure  5.  Profile  J  of  helium  density  from  three  different 
pair  potentials  surrounding  a  dopant  at  center  of  300 
atom  nanodroplct.  The  density  for  a  droplet  without  a 
dopant  is  also  shown. 

2.3  Solvation  of  A1  atoms  in  superfluid 
helittin  droplets 

2.3.1  A1  atom  excitation  spectra 

We  have  measured  the  LIF  spectrum  of 
the  3^D  <-  3^P  transition  of  A1  atoms  in 
superfluid  helium  nanodroplets  (see  Figure  6). 
The  spectrum  shows  the  A1  atoms  are  located 
inside  the  liquid  droplet,  having  a  large 
blueshift  (rou^ly  1000  cm  '  from  the  gas- 
phase  transition)  and  being  quite  broadened  (fit 
to  two  gaussians  with  FWHM  >400  cm'', 
separated  by  >  500  cm'').  The  observed 
splitting  is  much  larger  than  the  gas-phase 
3^P3/2  -  3^Pi/2  level  separation;  we  have 
performed  calculations  which  indicate  that  the 
orbital  angular  momentum  in  the  3^P  state  is 
quenched  by  the  helium  matrix  due  to  the 
strongly  anisotropic  interactions  with  the 
helium. 

23.2  A1  atom  time-resolved  measurements 

Lifetime  measurements  of  excited  A1 
atoms  (see  fig.  7)  were  accomplished  using  our 
time-correlated  single  photon  counting 
apparatus.  Temporal  evolution  of  the  total 
fluorescence  from  excitation  of  the  3^D  <—  3^P 
transition  yielded  lifetimes  of  6.4  ±0.1  ns,  in 
marked  disagreement  with  the  known  lifetimes 


Wavenumber  (cm‘^) 

Figure  6.  The  3^D  3^P  transitions  of  Al  atoms 

solvated  in  helium  nanodroplets.  The  fit  to  the  spectrum 
(solid  line)  is  a  summation  of  two  Gaussian  lineshapes 
(dotted  lines). 

of  the  3^D5/2  and  3^D3/2  states  (13.5  ns  and  13.3 
ns,  respectively)^'. 

Using  a  band-pass  filter  to  select  only 
the  3^D  ->  3^P  emission  resulted  in  no 
fluorescence  being  collected,  whereas  a  filter 
passing  only  fluorescence  from  the  4  S  — >  3  P 
transition  yielded  data  with  the  same  6.4  ns 
lifetime.  Exciting  directly  to  the  4^S  state  also 
gives  a  6.4  ns  lifetime.  It  is  apparent  that  the 
initially-populated  3^05/2,3/2  states  of  Al  are 
quenched  to  the  4^S  state  on  a  sub-nanosecond 
time  scale  due  to  interaction  with  the  He 
matrix.  This  entails  a  transfer  of  nearly 
7100  cm''  of  energy  into  the  surrounding 
helium.  This  energy  corresponds  to  the 
evaporation  of  over  1400  helium  atoms  from 
the  cluster  in  order  to  re-establish  the  pre¬ 
quenching  cluster  temperature  of  0.37  K, 
assuming  that  all  7100  cm''  of  energy  goes  into 
evaporation.  These  data  have  been  published 
in  ref  15. 

2.3.3  Al-He  potential  energy  surfaces 

The  (X)l^n  and  l^X  potential  energy 
surfaces  of  AlHe  have  been  calculated  using 
the  Hartree-Fock  Damped  Dispersion  (HFD) 
ansatz^^.  The  AEscf  component  was  calculated 
in  the  6-311-h-G**  basis  using  the  Gaussian 
98^^  program  suite,  correcting  at  each  point  for 
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Figure  7.  Dotted  line:  Temporally  dispersed 

fluorescence  of  A1  atoms  upon  excitations  near  the 
3^D  — >  3^P  transition  of  AL  The  rising  edge  is  masked 
by  scattered  light  Solid  line:  Temporally  dispersed 
fluorescence  of  A1  atoms  upon  excitations  near  the 
3^D  3^P  transition  of  Al,  detecting  only  frequencies 

near  the  4^S  ^  3^P  transition  frequency  region  of  AL 
The  inset  highlights  the  rising  edge  (solid  curve)  and  the 
best  fit  to  the  data  (dotted  curve). 

basis  set  superposition  error.  The  correlated 
(AEcorr)  component  makes  use  of  the  AlHe  long 
range  coefficients  Cn  calculated  using  the  AlAr 
Cn  coefficients  of  LePicard  et  the  Ar2  and 
He2  Cn  coefficients  of  Standard  and  Certain^"^, 
and  standard  combination  methods^^.  These  Cn 
coefficients  used  in  generating  the  1  11  and  1  Z 
AlHe  HFD  surfaces  are  listed  in  Table  6. 

Spin-orbit  coupling  of  the  l^H  and  l^Z 
AlHe  HFD  surfaces  yields  the  solid  curves 
given  in  Figure  8,  again  with  energies  referred 
to  the  original  Al([Ne]3s^3p*)  +  He(ls^) 
asymptote. 

When  spin-orbit  interaction  is  included, 
the  l^n  surface  is  found  to  split  into  a  l^ni/2 
and  1^113/2  surface,  the  former  of  which  mixes 
with  the  l^Zi/2  state  and  thus  has  a  minimum  of 


Figure  8.  Spin— orbit  averaged  1^11  and  and  spin- 
orbit-decoupled  l^Z|/2,  l^ni/2,  and  1^03/2  potential 

energy  surfaces  of  AlHe. 


parameter 

1^1  Al-He 

l^n  Al-He 

C6 

1.6695x10^ 

8.6183x10'* 

Cg 

2.1140x10^ 

1.0913x10^ 

Cio 

3.7354x10’ 

1.9283x10’ 

Table  6.  coefficients  used  in  generating  the  l^n  and 
l^Z  Al-He  HFD  surfaces.  All  energies  in  cm"'  and 
distances  in  A. 


potential  at  larger  intemuclear  distance 
(5.05  A)  than  either  the  or  the  spin-orbit 
averaged  1^11  surface.  Characteristic 
parameters  of  these  surfaces  are  given  in  Table 
7.  The  l^n  surface  is  seen  to  have  a  minimum 
at  4.12  A  while  the  l^Z  minimum  occurs  at  an 
AlHe  distance  of  5.77  A. 

After  generation  of  our  HFD  AlHe 
surfaces,  we  have  used  them  to  investigate  the 
nature  of  Al  electron  localization  for  the  Al- 
Hen  system.  In  particular,  the  main  question 
addressed  here  is  whether  the  ground  state  Al 
valence  is  best  described  by  the  spin-orbit 
decoupled  1^113/2,1/2  surfaces  or  by  the  spin- 
orbit  averaged  1^11  surface  (aspherical).  In  the 
former  case,  we  expect  only  the  lower-energy 
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Rmm(A) 

efcrn') 

ofcm'') 

i^n 

4.12 

16.45 

3.64 

l^Z 

5.77 

3.74 

5.1 

1^03/2 

4.12 

16.42 

3.64 

i^n,/2 

5.05 

4.97 

4.32 

l'll/2 

5.54 

3.81 

4.94 

Table  7.  Parameters  of  the  AlHe  potenetial  energy 
surfaces  as  calculated  by  the  HFD  ansatz. 


(l^ni/2)  surface  to  be  populated  due  to  the  cold 
environment  of  the  helium  droplet.  The  A1 
electron  density  correlating  with  this  state  is 
spherical.  If  the  spin-orbit  interaction  is 
quenched  in  the  helium,  however,  the  spin- 
orbit  averaged  l^n  surface,  with  an  aspheric  A1 
electron  density,  would  describe  the  Al-Hen 
interaction.  ITie  upper  frame  a  of  Figure  9 
compares  the  AlHen  (n  =  1,  4,  6,  7,  8)  1^11 
potential  energy  surfaces  based  upon  our 
1  ni/2AlHe  HFD  results  used  under  the 
assumption  of  additivity. 

These  surfaces  are  created  by  arranging 
the  helium  atoms  in  symmetric  polygons  in  the 
nodal  plane  aroimd  the  p  orbital  of  Al.  From 
the  figure  it  is  seen  that  a  minimum  appears  at 
distances  which  are  near  the  initial  spin-orbit 
averaged  1  11  minimum  (4.13  A),  while  the 
l^ni/2  minimum  at  5.05  A  goes  through  a 
“shelf’  state  and  almost  disappears.  We  can 
thus  conclude  that  as  the  number  of  He  atoms 
around  the  node  of  the  Al  atom  increases  the 
mixing  of  Z  character  into  the  l^H  state, 
predicted  by  spin-orbit  mixing,  becomes  less 
and  less.  This  lack  of  Z  character  corresponds 
to  less  Al  valence  electron  density  in  the  xy 
plane  (i.e.,  the  pleme  defined  by  the  Hen 
polygons),  creating  an  attractive  region  into 
which  the  He  atoms  can  be  drawn.  Another 
way  to  view  these  changes  in  the  potential 
energy  surfaces  is  as  representative  of  the 
localization  of  the  Al  valence  electron  in  its  pz 
orbital  (/.e.,  the  Al  p  orbital  perpendicular  to 
the  plane  defined  by  the  He  atoms),  which 
occurs  to  a  greater  extent  as  n  increases.  Thus 
these  additive  surfaces  support  the  conclusion 
that  the  spin-orbit  decoupling  is  quenched  in 
the  helium  solvent  so  that  the  density  profile  of 


Figure  9.  (upper  panel)  l^n  potential  energy  surfaces 
of  AlHe  (dot-dot-dash  line),  AlHe4  (dot-dash  line), 
AlHcfi  (dotted  line),  AlHe7  (d^ished  line),  and  AlHcj 
(solid  line),  where  the  helium  atoms  are  arranged 
symmetrically  on  the  xy  plane  around  the  Al  Pz  orbital, 
(lower  panel)  Probability  of  occupation  of  the  Al  Pz 
orbital  as  a  function  of  the  Al-Ke  distance  for  each  of  the 
AlHe„  configurations  depicted  above,  using  the  same 
line  designations. 

the  Al  p  electron  takes  on  the  aspherical  pz 
character. 

The  lower  frame  b  of  Figure  9  shows 
the  probability  of  occupation  of  the  pz  orbital 
as  a  function  of  Al-He  distance  for  &e  AlHon 
(n  =  4,  6,  7,  8)  surfaces  depicted  in  the  upper 
frame  of  the  Figure.  The  distance  coordinate 
corresponds  to  a  symmetric  “breathing”  of  the 
polygonal  He  ring.  It  is  found  in  comparing 
the  two  frames  that  the  movement  from  spin- 
orbit  averaged  (0.33)  population  of  the  Al  pz 
orbital  to  full  (1,0)  occupancy  occurs  in  the 
region  between  the  two  minima  and  at  a  longer 
intemuclear  distance  as  the  number  of  He 
atoms  in  the  node  is  increased.  It  is  thus  seen 
that  spin-orbit  mixing  is  more  rapidly 
quenched  as  the  munber  of  He  atoms  forming 
the  polygon  (/.«.,  as  the  number  of  He  atoms 
surrounding  the  Al  on  its  nodal  plane) 
increases.  The  He-He  distance  corresponding 
to  the  peak  of  the  low-temperature  pair 
correlation  of  liquid  helium  (3.5  A)^^  provides 
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the  upper  limit  for  the  He-He  distances  of  the 
He  atoms  which  make  up  the  Al-centered 
polygons^^.  The  Al-He  distance  of  4.2  A 
corresponds  both  to  the  AlHe  1^11  equilibrium 
intemuclear  distance  and  to  the  position  of  the 
AlHen  ground  state  minimum  for  /i>4  (cf. 
Figure  9).  At  this  AlHe  distance,  a  ring  of  7(8) 
He  atoms  is  found  to  have  a  He-He  spacing  of 
3.6(3.2)  A.  Thus  a  polygon  containing  7  He 
atoms  exhibits  He-He  distances  which  are  too 
large.  Therefore  we  conclude  that  at  least  8  He 
atoms  will  collect  around  the  A1  atom.  This 
has  been  shown  above  to  be  sufficient  to 
quench  the  spin-orbit  coupling. 

2»4  Theory  of  intermolecular  forces  and  its 
application  to  light-metal  atom  complexes. 

The  Hartree-Fock  plus  Dispersion 
method  has  two  components:  the  uncorrelated 
energy,  calculated  at  the  Hartree-Fock  self- 
consistent  level,  and  the  correlated  energy, 
calculated  by  an  inverse  distance  multipole 
expansion.  Each  term  in  the  expansion  is 
individually  damped  to  account  for  charge- 
overlap.^  These  damping  functions  are 
modeled  on  the  hydrogen  triplet  dimer,  which 
can  be  calculated  exactly.  The  size  of  the 
atoms  is  used  to  scale  the  amount  of  overlap  by 
taking  the  inverse  ratio  of  the  ionization 
potential  with  hydrogen  and  raising  it  to  0.66 
(so  p=(IP"/IP^)°®^),  where  p  is  the 
dimensionless  distance  scaling  parameter  and 
X  is  the  atom  or  molecule  of  interest).  This 
method  does  not  include  coupling  between  the 
exchange  and  correlation  energies. 

For  the  rare  gas  atoms,  the  HFD 
method  has  produced  good  potentials  for  the 
dimers^  (errors  in  binding  energy  and 
equilibrium  distance  of  1%  or  less).  Yet  it  is 
well-known  that  this  method  fails  when 
“softer”  atoms  (with  a  low  IP,  i.e.,  alkalis)  are 
involved.  Because  of  the  importance  of  “soft” 
atoms  to  the  HEDM  program  we  have  studied 
how  to  improve  the  predictive  ability  of  the 
HFD  method  for  these  atoms  and  their 
mixtures. 


2.4.1  Dimers  of  Soft  Atoms 

For  the  alkali  triplet  dimers,  the  HFD 
error  is  large  (30%  or  more)  and  increases 
down  the  group,  always  yielding  too  attractive 
potentials.  The  alkaline  earth  dimers  are  also 
calculated  to  be  too  attractive.  The  valence 
shells  of  these  dimers  have  significantly  more 
overlap  than  the  corresponding  hydrogen  case. 
The  extra  attraction  is  a  result  of  the  HFD 
method  treating  them  like  undistorted  spheres 
centered  on  each  nucleus. 

These  are  “soft”  atoms,  however,  with 
easily  distortable  electron  clouds.  As  the 
nuclei  get  closer  together,  the  electrons  resist 
overlapping  and  are  moved  off-center. 
Effectively,  the  nuclei  continue  to  approach 
each  other,  while  the  electron  clouds  remain 
farther  apart.  So,  in  the  standard  HFD  method, 
the  multipole  expansion  is  calculated  at  a 
distance  that  is  too  short,  producing  too 
attractive  results. 

As  the  electron  clouds  are  pushed  away 
from  their  respective  nuclei,  a  quadrupole 
develops.  Using  the  quadrupole  moment 
calculated  in  Gaussian  98^'  and,  treating  the 
system  as  four  point  charges,  we  can  determine 
the  distance  between  the  nuclei  and  the  centers 
of  the  electron  clouds.  The  alkali  triplet  dimers 
and  Mg2  show  a  “lag”  of  the  intemuclear 
distance  of  about  0.5  bohr.  We  refer  to  this 
additional  term  as  the  quadmpole  correction, 
and  the  new  method  for  calculating  potentials 
as  HFDQ.  This  is  a  first  step  towards 
including  exchange-correlation  coupling  terms. 

Using  the  new  HFDQ  method  we  have 
greatly  reduced  the  error  in  the  alkali  triplet 
(and  magnesium)  dimer  potential  curves  to 
15%  or  less  (in  the  well  depth;  to  1%  for  the 
distance).  The  remaining  discrepancy  is 
probably  due  to  the  shape  of  the  distorted 
electron  clouds;  HFDQ  still  treats  the  clouds  as 
spheres.  The  method  works  as  well  for 
mixtures  of  alkalis,  for  instance  for  the  NaLi 
triplet  (see  table  8).  When  applied  to  alkali 
oligomers  (e.g.,  trimers  and  tetramers)  the 
errors  decrease  to  about  10%  (see  table  9). 
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kMM 

Li2 

4.19 

ISI 

4.175 

HFD 

461.5 

3.88 

Naz 

lit. 

174.96 

5.09 

HFDQ 

199.84 

5.01 

4.7 

Kz 

mm 

5.77 

HFDQ 

269.489 

HFD 

445.92 

HE! 

ESS! 

429.6 

ItiiM 

HFDQ 

3.95 

HFD 

702.65 

3.66 

215.9 

4.76 

HFDQ 

252.974 

4.61 

HFD 

337.06 

4.28 

Table  8.  Alkali  triplet  dimer  potential  well  parameters 
from  the  HFD  and  HFDQ  methods,  compared  to 
literature  values. 

The  concept  of  distorted  electron 
clouds  also  accounts  for  the  non-additivity  of 
the  SCF  repulsion  Which  turns  out  to  be 
responsible  for  more  than  90%  of  the  additivity 
of  the  whole  Naa  potential  (see  the  end  of 
Section  2.1.2  and  ref.  12).  The  HF  energy  of 
the  trimer  and  tetramer  show  less  repulsion 
than  adding  the  SCF  of  the  dimer  pairs.  With 
the  softer  atoms,  the  nuclei  are  able  to 
approach  closer  because  the  electrons  are 
pushed  out  of  ftie  way. 

2.4.2  Mixtures  of  Soft  and  Hard  Atoms 

Neither  HFD  nor  HFDQ  work  for 
mixtures  of  soft  and  hard  atoms  (e.g.,  the 
alkali-rare  gas  dimers).  Most  of  the  potentials 
have  wells  too  shallow,  but  a  few  are  just  right, 
or  even  too  deep.  Why  does  the  method  fail? 
Could  it  be  the  difference  in  sizes  of  the 
atoms?  Not  likely,  since  LiH  is  too  deep,  while 
NaH  is  too  shallow;  LiHe  is  correct,  while 
NaHe  is  again  too  shallow.  And  both  alkali 
complexes  with  neon  are  much  too  shallow 
(see  table  10). 

Probably  size  is  important,  but  only 
through  the  distortion  of  the  electron  cloud.  As 
the  hard  partner  approaches,  the  soft  atom 
distorts  around  it,  almost  forming  a  dimple  in 


Dimer 
e/cm"'  Rffl/A 

Trimer 
e/cm"'  Rm/A 

Tetramer 
e/cm‘'  Kjk 

Na  lit. 

175.0 

5.09 

835.8 

4.40 

3.9 

HFDQ 

199.8 

5.01 

962.4 

3566 

3.83 

HFD 

KB 

5200 

3.64 

KB 

8600 

3.1 

IBI 

3.95 

KeJilK 

3.55 

8100 

3.15 

HFD 

3.66 

KKsI 

3.22 

14140 

2.95 

Table  9.  Comparisons  of  potential  well  parameters  for 
oligomers  of  Na  and  Mg. 


LiH 

e/cm'' 

RJA 

NaH 

e/cm"' 

RJA 

iiifcM 

HFDQ 

4.87 

5.99 

HFDQ 

3.95 

6.32 

5.98 

^1 

LiHe 

NaHe 

HFDQ 

1.50 

6.12 

HFDQ 

1.18 

6.54 

HFD 

1.51 

6.13 

HFD 

1.18 

6.52 

LiNe 

NaNe 

lit. 

8.78 

5.12 

lit. 

8.0 

5.3 

HFDQ 

5.84 

5.4 

HFDQ 

4.64 

5.70 

HFD 

6.00 

5.36 

HFD 

4.68 

5.68 

LiAr 

NaAr 

lit. 

42.2 

4.89 

lit. 

41.7 

5.01 

HFDQ 

34.7 

5.02 

HFDQ 

34.8 

5.05 

HFD 

42.4 

4.72 

HFD 

37.8 

4.86 

LiKr 

NaKr 

lit. 

68.0 

4.78 

lit. 

68.9 

4.74 

HFDQ 

45.3 

5.12 

HFDQ 

44.8 

5.22 

HFD 

57.8 

4.74 

HFD 

53.7 

4.85 

Table  10.  Comparisons  of  potential  well  paramters  for 
several  sofl/hard  atom  combinations. 


which  the  hard  atom  sits.  The  distortion  then 
slightly  polarizes  the  hard  partner  allowing  for 
increased  correlation  energies  and  greater 
attraction.  Judging  by  the  worse  results 
obtained  with  HFDQ  it  is  clear  the  shape  of  the 
distortion,  and  not  just  the  distance,  is 
important. 

Further  efforts  at  determining  the  shape 
of  the  electron  cloud  will  improve  our 
understanding  of  the  alkali-alkali  interaction, 
the  alkali-noble  gas  mixtures  and,  ultimately, 
of  the  alkali-hydrogen  interactions. 
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2.5  Progress  towards  assembline  a  facility 
for  stimulated  Raman  spectroscopy  (SRS)  in 
doped  and  undoped  hydrogen  clusters. 

2.5.1  The  SRS  experiment 

In  recent  years,  the  experimental 
feasibility  of  observing  a  superfluid  phase  of 
para-hydrogen  (P-H2)  has  been  investigated  in 
various  contexts;’*’^’’^®'^'’^^  to  date,  the  existence 
of  such  a  phase  has  not  been  conclusively 
proven. 

The  goal  of  an  ongoing  effort  here  at 
Princeton  is  to  produce  ultracold  p-H2  droplets 
of  a  few  hundred  to  a  few  thousand  molecules 
in  a  supersonic  expansion,'*^  which  will  be 
spectroscopically  investigated  using  SRS.  Such 
droplets  have  good  prospects  for  being 
superfluid  because  both  their  small  size  and 
their  very  short  lifetime  in  the  experiment 
decrease  the  chance  of  freezing  nucleation:  it 
may  then  be  possible  for  them  to  be 
supercooled  so  much  below  the  triple  point  of 
P-H2  that  Bose-Einstein  condensation  occurs. 

Spectroscopic  evaluation  of  the 
possibly  superfluid  aggregation  state  of  the 
hydrogen  nanodroplets  calls  for  Raman 
scattering,  since  infrared  transitions  in 
condensed  molecular  hydrogen  are  very  weak, 
induced  only  by  lattice  defects,  surfaces,  and 
impurities.  Due  to  the  low  number  density  of 
hydrogen  molecules  in  a  cluster  beam, 
stimulated  Raman  scattering  is  called  for 
because  of  its  increased  eflHciency.  During 
stimulated  Raman  scattering,  a  high-energy 
photon  (738  nm  in  our  case)  is  absorbed  while 
at  the  same  time  emission  of  a  lower-energy 
photon  (1064  nm)  is  stimulated.  The  difference 
in  energy  promotes  a  hydrogen  molecule  to  the 
v=l  vibrationally  excited  state.  Subsequent 
relaxation  back  to  the  ground  state  leads  to 
partial  evaporation  of  the  hydrogen 
nanodroplet.  As  in  the  IR  spectroscopy  of 
doped  ^He  droplets,'*^  the  decrease  in  beam  flux 
can  be  detected  with  a  cold  bolometer. 
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Figure  10.  Spectral  regions  accessible  with  the  Indigo 
TiiAbOs  in  combination  with  the  Antares  Nd:YAG.  The 
label  (m,n)  refers  to  mixing  the  m'*'  harmonic  of  the 
Nd;YAG  with  the  n*  harmonic  of  the  Ti;Al203.  For 
example,  (-1,3)  is  the  difference  frequency  between  the 
tripled  Ti:Al203  and  the  fundamental  of  the  Nd:YAG. 
The  bottom  line  summarizes  the  accessible  wavelength 
ranges. 

2.5.2  The  laser  setup 

The  laser  system  for  stimulated  Raman 
scattering  in  hydrogen  consists  of  a 
continuously  tunable  Ti:Al203  “Indigo”  laser 
that  produces  narrow  band  radiation  around 
738ran,  and  a  Q-switched,  cavity-dumped 
Nd:YAG  laser  “Antares”  operating  at  1064nm. 
The  short  pulses  produced  by  the  Indigo  can 
easily  be  frequency  doubled  and  tripled;  when 
combined  with  the  Antares  and  its  harmonics, 
a  very  large  range  of  frequencies  can  be 
covert,  in  which  a  wealth  of  high-resolution 
experiments  can  be  performed  (see  Figure  10). 

2.5.3  The  Indigo  Ti;Al->0^  cavity 

The  narrow-band  output  of  the  Indigo 
is  produced  in  a  modified  Fox-Smith 
resonator'**  consisting  of  a  high  reflector  (HR), 
an  output  coupler  (OC),  a  grating,  and  a  tuning 
mirror  (see  Figure  11).  There  are  two 
interfering  light  paths  in  such  a  resonator: 
between  the  HR  and  the  OC  via  specular 
reflection  off  the  grating  (path  A),  and  between 
the  HR  and  the  tuning  mirror  via  first-order 
diffraction  off  the  grating  (path  B).  The  latter 
path  is  wavelength  selective,  because  the 
combination  of  grating  and  tuning  mirror  only 
reflects  rays  of  a  very  narrow  wavelength 
range  back  onto  themselves. 
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738  nm  (tunable)  +  1064  nm  for  SRS  in  hydrogen 


^  harmonic 
£  generation 


245  -  275  nm  tunable 


Figure  11.  Setup  of  the  laser  system.  The  Indigo  is  a 
tunable,  pulsed  single-mode  Ti:Al203  laser  that  can  be 
frequency-doubled  and  -tripled.  The  Antares  is  a  pulsed 
Nd:YAG  laser  synchronized  with  the  Indigo,  to  be  used 
for  stimulated  Raman  scattering  or  frequency  mixing 
with  the  Indigo.  The  diagnostics  are  used  both  for 
wavelength  measurements  and  mode  stabilization  of  the 
Indigo, 

The  lasing  threshold  of  the  Indigo 
resonator  is  greatly  reduced  if  we  have 
constructive  interference  between  the  two  light 
paths.  If  all  other  cavity  elements  are  fixed,  we 
can  change  the  length  of  path  A  by  displacing 
the  OC,  which  is  done  by  applying  a  voltage  to 
the  OC  PZT.  This  shifts  the  absolute 
wavelengths  at  which  constructive  interference 
occurs,  and  we  can  thus  make  sure  that  at  the 
maximvim  of  the  wavelength  selectivity  of  path 
B  we  simultaneously  have  constructive 
interference  of  the  two  light  paths. 

Furthermore,  we  must  place  the  HR  at  a 
position  where  both  light  paths  simultaneously 
have  a  node,  so  that  the  two  light  paths  actually 
form  an  eigenmode  of  the  laser  cavity.  We  do 
this  by  applying  a  voltage  to  the  HR  PZT.  This 
is  a  modification  to  the  original  design  of  the 


Indigo  by  Positive  Light,  Inc.,  in  which  there 
was  no  provision  for  such  fine-tuning.  We 
have  foimd,  however,  that  single-mode 
operation  of  the  Indigo  cannot  be  achieved 
without  this  additional  degree  of  freedom. 

We  have  put  in  place  a  computer-based 
feedback  system  that  observes  the  spectral 
quality  of  each  laser  shot  of  the  Indigo  through 
a  plane-parallel  6talon  (FSR  4.4  GHz)  imaged 
with  a  linear  camera,  and  then  adjusts  the 
voltages  on  the  two  cavity  PZTs  accordingly. 
Misplacement  of  the  HR  is  detected  by 
observing  the  emergence  of  neighboring  cavity 
modes,  spaced  by  about  800  MHz  from  the 
main  laser  output  frequency.  Misplacement  of 
the  OC  is  detected  by  the  emergence  of 
neighboring  maxima  of  constructive 
interference,  spaced  by  about  3  GHz  from  the 
main  laser  frequency.  The  different  spacings  of 
cavity  modes  and  interference  maxima  easily 
allow  their  identification  in  the  feedback 
algorithm.  With  this  system  in  place,  single¬ 
mode  operation  of  the  Indigo  can  usually  be 
established  within  a  few  hundred  milliseconds. 
As  the  tuning  mirror  is  rotated  slowly  and  thus 
the  wavelength  selected  by  path  B  is  scanned, 
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the  positions  of  the  HR  and  OC  mirrors  are 
adjusted  continuously,  allowing  smooth  tuning 
of  the  single-mode  output  of  the  Indigo.  The 
linear  tunability  is  now  only  limited  by  the 
gain  medium. 

After  stabilization,  the  Indigo  puts  out 
single-mode  radiation  that  is  essentially 
transfer  limited  by  the  pulse  duration  of 
approximately  19  ns.  Pumped  by  a  10  W 
“Evolution  X”  diode-pumped,  intra-cavity 
fi'equency-doubled  Nd:YLF  laser,  the  Indigo 
can  deliver  up  to  1.5  mJ  pulses  around  800  nm 
at  1  kHz  repetition  rate. 

The  wavelength  of  the  Indigo  radiation 
is  monitored  ‘vith  a  Burleigh  WA-4500 
wavemeter.  For  absolute  calibration,  we  have 
furthermore  measured  the  two-photon 
absorption  spectrum  of  rubidium  vapor  around 
778.1  nm  (see  Figure  12).  Combined,  these 
measurements  determine  the  Indigo  frequency 
to  within  less  than  100  MHz. 

2.5.4  Indiao  Harmonics  Generation 

The  higii  peak  power  of  the  Indigo 
invites  generation  of  second-  and  third- 
harmonic  laser  radiation.  Second-harmonic 
generation  (SHG)  is  done  by  focusing  the 
fundamental  into  an  angle-tuned  LBO  crystal, 
and  produces  up  to  200  |iJ  pulses  around  400 
nm.  Third-harmonic  radiation  is  produced 
(THG)  by  co-focusing  the  fundamental  and  the 
second  harmonic  into  an  angle-tuned  BBO 
crystal.  Since  the  polarization  of  the 
secondharmonic  is  orthogonal  to  the 
fundamental  due  to  bireMngence  phase 
matching  in  the  LBO  crystal,  we  rotate  the 
plane  of  polarization  of  the  fundamental  by  90° 
with  a  half-wave  plate  (HWP)  after  separating 
it  fi'om  the  second  harmonic  with  a  dichroic 
mirror.  We  can  produce  up  to  3.5  pJ  of  UV 
light  at  260  nm  in  pulses  of  about  7  ns  length. 

We  expect  to  receive  the  components 
for  fourth-harmonic  generation  (FHG)  soon. 
This  will  allow  us  to  extend  the  range  of 
accessible  wavelengths  all  the  way  down  to  the 
limit  of  vacuum  UV. 

The  angles  of  all  the  harmonic- 
generation  crystals  can  be  adjusted  by 


Figure  12.  The  two-photon  fluorescence  spectrum  of 
gas-phase  rubidium,  as  measured  with  the  Indigo.  The 
main  peak  corresponds  to  the  5Si/2  (F=3)  5D5/2  (F-5) 


transition  in  Rb,  and  is  at  385  285  142  378.280  ±  2 
kHz.‘“ 


picomotor  screws,  which  allow  software 
control  of  the  phase  matching  condition  as  the 
Indigo  is  scanned  in  frequency. 

2.5.5  Antares  Nd;YAG  Cavity 

We  have  converted  a  Coherent  Antares 
cw  Nd:YAG  laser  to  pulsed  operation,  in 
which  it  can  produce  10  ns  pulses  at  a 
repetition  rate  of  up  to  3  kHz,  and  with  a  pulse 
energy  of  up  to  2  mJ  (measured  at  1  kHz).  This 
was  achieved  by  inserting  into  the  cavity  a 
KD*P  Pockels  cell  (PC),  a  quarter-wave  plate 
(QWP),  and  a  thin  film  polarizer  (TFP),  in  the 
setup  schematized  in  Figure  10.  Switching  the 
Pockels  cell  on  starts  the  buildup  of  laser 
radiation  in  the  cavity,  and  switching  it  off 
dumps  the  entire  laser  power  out  of  the  cavity 
within  about  one  round-trip^  time.  This 
“dumped”  pulse  can  be  synchronized  with  the 
Indigo  pulse  to  within  less  than  2  ns,  which  is 
sufficient  for  stimulated  Raman  scattering 
experiments. 

The  unseeded  Q-switched  Antares  has  a 
spectral  width  of  about  0.6  cm"',  which  is 
unacceptable  for  high-resolution  Ramm 
spectroscopy.  To  remedy  this,  we  have  added  a 
Lightwave  S-100  seed  laser  that  reduces  the 
line  width  of  the  Antares  to  single-  or  double¬ 
mode  operation.  A  circuit  for  locking  the 
Antares  cavity  length  to  the  seeder  wavelength, 
to  produce  single-mode  output,  remains  to  be 
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implemented.  However,  the  cavity  modes  of 
the  Antares  are  spaced  by  only  140  MHz,  and 
therefore  operation  on  two  modes  introduces  a 
frequency  imcertainty  comparable  to  the  Indigo 
line  width,  and  further  narrowing  may  turn  out 
to  be  unnecessary. 
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01.2  Farady  Discussions  118:  Cluster  Dynamics,  Durham,  U.K.,  April  18-20,2001. 

01 .3  APS  Division  of  Atomic,  Molecular,  and  Optical  Physics  Meeting,  London,  Ontario, 

May  16-18,2001. 

01.4  Air  Force  Office  of  Scientific  Research  Contractors  Meeting,  Irvine,  CA.,  May  21-23, 2001 

01 .5  Fifty-sixth  International  Symposium  on  Molecular  Spectroscopy,  Columbus,  OH., 

June  11-15, 2001. 


5.3  Invited  seminars  given  by  G.  Scoles  in  which  work  carried  out  under  this  grant  was 
discussed.  frotal  =  24) 


99.1  Phys.  Dept.,  University  of  Delaware,  Feb.  24 

99.2  University  of  Florence,  Italy,  March  16 

99.3  University  of  Bologna,  Italy,  March  1 

99.4  University  of  Genova,  Italy,  March  3 

99.5  Stevens  Institute  of  Technology,  Hoboken,  N.  J.,  Nov.  17 

99.6  University  of  Wisconsin,  Chemistry,  Nov.  30 
00. 1  Phys.  Dept.  SISSA,  Trieste  (Italy)  Jan.  1 1 
00.2  Chem.  Depi.,  Ohio  State  University,  Feb.  23 
00.3  Chem.  Dept.,  Ohio  State  University,  Feb.  25 
00.4  Phys.  Dept ;  Harvard,  March  8 

00.5  Phys.  Dept.,  U.of  Barcelona  (Spain),  April  26 

00.6  Nat.Res.Council  of  Canada,  Ottawa,  May  12 

00.7  Phys.  Dept.,  U.  of  Genova  (Italy),  June  14 

00.8  Chem.Dept.  Tech.  Univ  of  Munich  (Germany),  June  30 

00.9  Chem.Dept.,  Utah  State  University,  Oct.  23 

01.1  ETH  Zurich,  Dept,  of  Chemistry,  May  9 

01 .2  Dept,  of  Chem.,  Univ.  of  Geneva,  May  1 5 

01.3  Dept,  of  Phys.Chem.,  Univ.  of  Basel,  May  1 6 

01 .4  Dept,  of  Phys.  Chem.,  EPFL,  Lausanne,  May  1 7 

01.5  Dept,  of  Phys.  Chem.,  EPFL,  Lausanne,  May  18 

01.6  Dept,  of  Phys.Chem.,  Univ.  of  Bern,  May  22 

01.7  Dept.  Phys.  Chem.,  Univ.  of  Helsinki,  June  18 

01.8  Dept,  of  Chem.,  Univ.  of  Alberta,  Edmonton  (Canada)  Nov.21, 22 


01 .9  Dept,  of  Physics,  Perm  State  Univ.,  College  Park,  PA,  Dec.  7 


5.4  Invited  seminars  given  by  Keviii  Lehmann  in  which  work  carried  out  under  this  grant 
was  discussed.  HTotal  =  8) 

97. 1  Department  of  Chemistry,  New  York  University,  Dec.  19, 1 997. 

98.1  Department  of  Chemistry,  Georgia  Technical  University,  Atlanta,  GA.,  Oct.  29, 1998. 

98.2  Department  of  Chemistry,  Emory  University,  Atlanta,  GA.,  Oct.30,  1998. 

98.3  Department  of  Physics,  University  of  Rhode  Island,  Dec.  11, 1998. 

99.1  George  Harrison  Spectroscopy  Laboratory,  MIT,  Oct.  19, 1999. 

01 . 1  Department  of  Chemistry,  Harvard  University,  Feb.  1 5, 2001 . 

0 1 .2  Department  of  Chemistry,  University  of  Southern  California,  Feb.  26, 200 1 . 

01 .3  Department  of  Chemistry,  University  of  Wisconsin  at  Eau  Claire,  Oct.  19,  2001 

5.5  Contributed  papers  at  national  and  international  meetings  by  other  members  of  the  £roup 
(Total  ==12) 

J.  Higgins:  Helium  Cluster  Isolation  Spectroscopy  of  Alkali  Dimers  in  the  Triplet  Manifold.  Gordon 
Conference  on  Molecular  and  Ionic  Clusters.  Ventura,  CA.  Jan.  4-9,  1998.  (poster  presentation) 

J.  Reho:  Time-resolved  Spectroscopy  ofNa.  Atoms  and  Oligomers  on  the  Surface  o/He  Clusters. 

The  13*  Symposium  on  Chemical  Physics,  Waterloo,  ON,  Canada.  Oct.  24-26,  1997  (poster 
presentation) 

rH 

J.  Reho:  Non-adiabatic  Processes  in  the  Photodissociation  of  Quartet  State  Naa  and  K3.  53  Ohio 
State  University  International  Symposium  on  Molecular  Spectroscopy.  Columbus,  OH.  June 
15-19,1998.  (oral  presentation) 

J.  Reho:  Direct  State-selective  Measurement  of  Cooling  Rates  of  Vibrationally  Excited  Na^ 

on  the  Surface  of  Helium  Clusters.  53"*  Ohio  State  University  International  Symposium  on 
Molecular  Spectroscopy.  Columbus,  OH.  June  15-19, 1998.  (oral  presentation) 

J.  Reho:  Dynamics  of  Optically-excited  Alkali  Metal  Dimers  and  Trimers  on  the  Surface  of  Helium 
Clusters.  216*  ACS  National  Meeting.  Boston,  MA.  Aug.  23-27,  1998.  (oral  presentation) 

M.  Radcliff:  Charge  Distribution  Deformation  in  van  der  Waals  Complexes  Involving  Atoms  With 
Low  Ionization  Potentials.  17*  Symposium  on  Chemical  Physics,  Waterloo,  ON,  Canada.  Nov.  9- 
1 1 ,  2001 .  (oral  presentation) 
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M.  Radcliff:  Charge  Distribution  Deformation  in  van  der  Waals  Complexes  Involving  Atoms  With 
Low  Ionization  Potentials.  2"**  PUUDel  Workshop  on  Intermolecular  Interactions,  Newark,  DE. 
Oct.  12, 2001.  (oral  presentation) 

M.  Radcliff;  Spectroscopy  o/Mg  Atom-Doped  Helium  Nanodroplets.  55*’’  Ohio  State  University 
International  Symposium  on  Molecular  Spectroscopy.  Columbus,  OH.  June  12-16, 2000.  (oral 
presentation) 

M.  Radcliff:  Spectroscopy  of  A1  Atoms  Solvated  in  Helium  Nanodroplets.  55*  Ohio  State 
University  International  Symposium  on  Molecular  Spectroscopy.  Columbus,  OH.  June  12-16, 
2000.  (oral  presentation) 

M.  Radcliff:  Spectroscopy  of  Atoms  in  Helium  Nanodroplets.  U*  PUUDel  Workshop  on 
Intermolecular  Interactions,  Princeton,  NJ.  June  15, 2000.  (oral  presentation) 

M.  Radcliff:  Studies  o/Mg  Atoms  Solvated  in  Helium  Nanodroplets.  15*  Symposium  on  Chemical 
Physics,  Waterloo,  ON,  Canada.  Nov.  5-7, 1999.  (poster  presentation) 

M.  Radcliff:  Non-Adiabatic  Processes  in  the  Photodissociation  of  Quartet  State  Sodium  Trimer. 
13*  Symposium  on  Chemical  Physics,  Waterloo,  ON,  Canada.  Oct.  24-26,  1997.  (poster 
presentation) 


6«0  Honors  and  Awards 


During  the  period  of  tenure  of  this  grant  the  PI  has  received  the  following  honors  and  awards; 
1997  -  Elected  Fellow  of  the  Royal  Society  (UK) 

1999  -  SamuelMcElvainLecturer,  University  of  Wisconsin,  Madison 

2000  -  Honorary  Science  Doctorate  from  the  University  of  Waterloo  (Canada) 

2000  -  Elected  Foreign  Member  ofthe  Royal  Academy  of  Arts  and  Sciences  of  The 
Netherlands 

2000  -  The  2000  Chemical  Physics  Symposium  ofthe  University  of  Waterloo  (Canada) 

was  dedicated  to  the  P.I.  65*  Birthday. 

2001  Gunning  Lecturer  at  the  University  of  Alberta  (Canada) 

2002  Peter  Debye  Award  for  Physical  Chemistry  of  the  ACS 
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